Sigma meson in QCD sum rules using a two quark current with derivatives 
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We study the a meson in QCD sum rules using a two quark interpolating field with derivatives. 
In the constituent quark model, the a meson is composed of a quark and an antiquark in the relative 
p-wave state and is thus expected to have a larger overlap with an interpolating field that measures 
the derivative of the relative quark wave function. While the sum rule with a current without 
derivative gives a pole mass of around 1 GeV, the present sum rule with derivative current gives a 
mass of around 550 MeV and a width of 400 MeV, that could be identified with the a meson. 



o. 
o- 

(N ; 

u '■ 

Or 

<: 



^3- 
Or 

i ' 

<l> : 

43. 



> 

(N 

00 

o 

O 
oo 
o 



X 



PACS numbers: 12.38.Lg,11.55.Hx,14.40.Cs 

The existence of the light scalar meson or <t(600) had 
been controversial for a long time despite of its impor- 
tant role in chiral dynamics in QCD physics [l|, |2j]. It is 
only recently that the existence has been confirmed [j, [j] , 
through the careful reanalyses of the 7T-7T scattering phase 
shift 5] , and the findings of the a pole in the heavy parti- 
cle decays such asD-> 7T7T7r and T(35) — > T7T7T [y, 0,[a]. 

On the other hand^the physical content of the a meson 
is still controversial [9J. In addition to the usual picture 
of a qq state, there are several plausible candidates like a 
four quark state qqqq [10, LUJ and a tttt molecule, which 
can be mixed with a glueball. Also, the a meson may be 
a collective oo states just as the it is, due to the chiral 
symmetry |12| . Hence, it is important to confirm the 
existence of the a meson, and to investigate its nature 
on the basis of QCD. 

In lattice QCD [13], most of the calculations using a 
two quark interpolating field seem to predict the ground 
state mass of scalar particle in the isospin 1 channel to 
be above 1.3 GeV 0, M, M^JM, while some pre- 
dict it to be around 1 GeV [19|, [20]. In contrast, most 
lattice calculations based on a four quark interpolating 
field consistently predict the mass to be around 1 GeV 
for the /o, a [lE|2l| and around 600 MeV for the a [H. 
The existence of the a with a two quark interpolating 
field was confirmed in a full QCD simulation for the first 
time in Ref. [23|. It was argued that the disconnected 
diagram is crucially responsible for reducing the a me- 
son mass to that comparable to the p meson mass for the 
lightest current quark mass achieved in the simulation. 
It is expected that if extrapolated to the chiral limit, the 
a meson mass would be such as m w < m a < m p . On the 
other hand, for the k, where no disconnected diagram 
contributes when a two quark interpolating field is used, 
the calculated mass was about two times that of the K* 
|22j ]. These results suggest that the a might have a very 
different quark structure than the rest of the scalar nonet. 

In this paper, we will study the a meson in QCD 



sum rules, which is another first principle approach to 
QCD [241. The first QCD sum rule attempt was made 
in Ref. [25] with the usual two quark interpolating field 
~ uu(x) 4- dd(x); a local operator. Assuming the de- 
cay width to be zero, the calculated mass was found 
to be degenerate with its isospin 1 partner and around 
m a = 1.00 ± 0.03 GeV, which is too large to be identify 
as the a meson. On the other hand, recent works using 
four quark interpolating fields [2y, [23, US, (2E [30] seem 
to give a mass closer to that of the physical a. 

Our idea in this work is to observe that in a non- 
relativistic quark model, the quark and anti-quark in the 
a meson are in the relative p-wave state, as can be de- 
duced from its quantum number J PC = ++ . The rel- 
ative wave-function has a node and a finite slope at the 
origin. Therefore, it is expected that the qq operator 
with additional covariant derivative would have a larger 
overlap with the physical a state [3l|, [32|, [33| • Thus we 
employ a two quark interpolating with derivatives, which 
will be called the non-local operator, 
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lu (\z ■ *d\ u{x) + d(iz- *D J d(x) } . 

(1) 

Here, \D = \D + (iDy, D^ = <9 M — \gA^ is the covariant 
derivative and z^ is a four- vector such that z 1 = — 1 in 
order to probe the wave-function in the space-like direc- 
tion. We note that for n = 0, it corresponds to the local 
operator. Of course, the quark operator in the interpolat- 
ing field is that of a relativistic current quark. Therefore, 
the current can still couple to the scalar meson that in 
the non relativistic limit is a p-wave state. Neverthe- 
less, we hope to increase the overlap with the physical 
a by employing an interpolating field that has the ex- 
pected non relativistic limit. Such interpolating currents 
with derivatives were previously used successfully to in- 
vestigate the p-wave nucleon excited states [32J, 33] . This 
choice of interpolating field would also be a natural choice 
in the lattice calculation, where the quark masses are still 
much larger than the physical limit. 

Let us start by considering the following time-ordered 
two-point correlation function 



|TO„(ar)O (0)|0). (2) 



T n0 (q ,z-q, z ) = i / d 



The correlation function will couple to the scalar meson 
when n is even. If n is too large, the correlation will be 
dominated by the continuum contribution. Therefore, 
the optimal choice for our analysis is n = 2. For n — 
2, the operator product expansion up to dimension 6 is 
calculated to be 
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Note that we have calculated the OPE up to dimension 6 
operator. It is important to calculate the OPE up to the 
so called quark tree diagram contributions when calculat- 
ing the OPE of current-current correlation function. This 
is so because purely gluonic operators are suppressed by 
loop effect that are typically suppressed by 1/47T, and 
higher order quark operators are suppressed by the ad- 
ditional factors of a s . This means that when calculating 
meson correlation function with multiple quark currents, 
the OPE should be investigated up to quark operators 
with quark number equal to 2 times that appearing in 
the current. In our case, since our current involves two 
quark operators, the minimal quark operator that has to 
be considered is the four quark operator of dimension 6. 
Picking up terms proportional to z 2 , subtracting out 
the continuum contribution that start from the threshold 
so, and then performing the Borel transformation, leads 
to the following expression. 

n (Af 2 ) = --Lm 6 E 2 ( Sq /M 2 ) - ^M 2 E ( So /M 2 )G Q 
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where M is the Borel mass, and 
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The theoretical side or the total OPE[24j, is defined as 
Eq.(j4| with sq —> oo. The sum rule for the a is ob- 
tained by equating Eq.(H]) to the spectral density below 
the threshold sq, which is dominated by the a, through 
the Borel transformed dispersion relation. 



n (M 2 ) = -/ dse- s/M2 Indies). 
n Jo 



(6) 



For the spectral density coming from the <r, we include 
a finite width by employing a simple Breit-Wigner form. 



Namely, 
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Here, / is the overlap constant of the current with the 
physical a. The right hand side of Eq.© with Eq.® 
will be denoted as 11% (M 2 ). 

The QCD sum rule for the width and mass of the a 
is obtained by taking the ratio of the Eq.© with its 
derivative with respect to the Borel mass. 
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where 
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. The left-hand side of Eq.® 



becomes to 2 in the narrow width approximation. In this 
limit, m 2 is first plotted as a function of M 2 for dif- 
ferent values of the threshold sq. From the plot with a 
threshold that gives the most stable plateau, the mass 
is determined from the value at the stable plateau re- 
gion. We remind the reader that such analysis has to be 
taken with the usual cautions as in any QCD sum rules 
involving ratios of derivatives 34 1 . 

In the present analysis, we will determine not only the 
mass but also the width, with the following prescription: 



1. Fix a continuum threshold sq, typically 0.< 
GeV 2 . 



2.0 



2. For each sq, calculate the Borel window M^ in and 
A^max by following the criterion that the power cor- 
rections and continuum correction should be re- 
spectively smaller then the theoretical side[24[ by 
certain fraction chosen as, 



M 2 



AP 



Power correction 

Total OPE sum 
Pole contribution 

Total OPE sum 



< 0.3-0.4 



> 0.3-0.5 



(10) 
(11) 



3. Compute the mass m a by solving Eq. ((5|) numeri- 
cally for each fixed T, which are varied from to 
500 MeV. 

4. Plot m a as a function of M 2 for each T, and de- 
termine to ct and r that gives the most stable curve 
within the Borel window. 

Now let us discuss numerical results. We begin with 
the n = case corresponding to the local operator [25j |. 
In Fig. [TJ we plot m a as a function of M 2 for various 
widths with sq = 1.6 GeV 2 . We see that with T = 0, we 



can reproduce the result of Ref . [25| , which shows a good 
Borel stability and gives the mass of about 1 GeV. How- 
ever, if we allow for a finite width, the stability is lost and 
we can not conclude that there is a resonance with that 
width. Therefore, it is quite unlikely that the physical 
a meson has strong coupling to the correlation function. 
Instead, the stable Borel curve with small width seems to 
suggest a stronger coupling to /o(980) or a (980), which 
are degenerate in the present calculation. 
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FIG. 1: Local operator: m a as a function of M 2 for various 
widths with so = 1.6 GeV 2 . Vertical dotted lines denote the 
Borel masses at which continuum contribution becomes 30%- 
50% and power contribution becomes 30%-40%, respectively. 

For the non-local operator case with n = 2, numerical 
analyses show that we can find a stable plateau with a 
maximum in the Borel curves within the Borel window 
for a wide range of the threshold parameter sq. The ob- 
tained mass and width are m a ~ 350 — 950 MeV and 
T - 400 - 500 MeV for s = 1.2 - 2.0 GeV 2 . In Fig. [2 
we show the result for so = 1.4 GeV 2 . We find that the 
Borel window is most stable and wide with the maximum 
within the window for sq = 1.4 ± 0.1 GeV 2 . It is very 
important to have a stable region or the extremum point 
within the Borel window. Otherwise, the sum rule loses 
any predictive power and suggests that it is dominated 
by the continuum. Even when s = 1.4 GeV 2 , we can see 
that such a stable Borel curve with a maximum appears 
when we introduce a large width. This is in contrast to 
the local operator case, and suggest that the the correla- 
tion function is saturated with a resonance with a large 
width. We note that the maximum Borel mass is deter- 
mined from the condition that the continuum contribu- 
tion is less than 70 %. Also, in this case, contribution of 
the highest dimensional power correction (|(g<?)| 2 ) is al- 
ways less than 10% within the Borel mass region shown 
here. 

Usually, meson sum rules are insensitive to the width. 
An example being the vector meson sum rule; while the 
OPE of the p and ui sum rule are identical their widths 
are very different. On the other hand, our sum rules are 
of higher dimensions than the vector meson sum rule. 
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FIG. 2: Non-local operator: m a as a function of M for var- 
ious widths with so = 1.4 GeV 2 . 



This means that the continuum contributes with higher 
moments and that the sum rules are more sensitive to 
the continuum values. This also means that the sum 
rules are more sensitive to the width. If the sum rule 
were of too high dimension, then the sum rule would be 
dominated by the continuum and we could not obtain any 
information about the pole or the width. What we have 
shown is that for the sum rules at hand, the dimensions 
are optimal such that the sum rules are sensitive to the 
width but at the same time still dominated by the pole, 
especially for the sum rule with non-local current. 

To further check the consistency in our approach, it 
would be useful to construct the three-point function 
among the two-quark current with derivative and the two 
pseudo-scalar currents, which will couple to the sigma 
and the pions, respectively. The analysis of such a func- 
tion will provide information for the a-ir-ir coupling, from 
which one could directly calculate the sigma width and 
check the consistency with the present calculation. This 
will be left as a future work. Another related work would 
be to introduce an effective diquark scalar field in QCD 
[40j | , construct the current for the sigma meson using the 
diquark current, and then perform a similar QCD sum 
rule analysis. The lowest dimensional current will not in- 
clude a derivative as two quarks and two anti-quarks in 
s-wave states will form a scalar meson in the constituent 
quark model. The QCD sum rule analysis will be more 
effective in concentrating on the attraction between the 
scalar diquarks and thus more effective than the previous 
QCD sum rule analysis using the four quark interpolating 
field. 

To summarize, we have investigated the a meson in 
QCD sum rules by using the non-local qq operator with 
a covariant derivative. Such a choice of current was mo- 
tivated by the fact that the wave-function is a p-wave 
in a naive non-relativistic quark model. We have con- 
firmed that with such a choice of current, the sum rule 
is consistent with a strong coupling to the a meson. We 
note that it is necessary to take into account both the 



mass and the width simultaneously: If we set the width 
to zero, we do not obtain the Borel stability at all. The 
best estimate that gives the most stable plateau is with 
a mass m a = 550 MeV and Y a = 400 MeV in the chiral 
limit, which is the first theoretical prediction based on 
a first principle QCD calculation. Our work also con- 
firms that the a has a large qq component. Therefore, 
it would be interesting to try similar approach in lattice 
calculation. There are several things to improve in the 
future. First, in the present analysis, the 1 = channel 
(a) and the 1=1 channel (S) are degenerate. To lift 
the degeneracy, we have to introduce higher dimensional 
operators or effectively take into account such contribu- 
tions through the instanton effect [35|, big, [37|. This is 
so because single instanton configurations are known to 
contribute to correlation function of scalar currents [38j]. 
Moreover, ao(980) is expected to have a large tetraquark 
component as has been shown in recent QCD sum rule 
calculations [2a, [23, [28|, [29J,[3(|. Secondly, it is necessary 



to take into account the a s correction to the perturbative 
contribution and the corrections from the finite current 
quark masses. Through such extension, we can inves- 
tigate the whole scalar nonet, among which the k(800) 
meson is of great interest [8|, |39| . Third, it would be inter- 
esting to study the other p-wave mesons s yste matically 
by using the non-local operator, as in Ref. [25j. 



Acknowledgments 

We would like to thank APCTP for sponsoring the 
workshop on 'Hadron Physics at RFHC' and the orga- 
nizers of NFQCD2008 at YITP. The discussions during 
these workshops have lead to the final version of this pa- 
per. We are particularly grateful to Teiji Kunihiro for 
useful discussions. 



[1] E. van Beveren, T. A. Rijken, K. Metzger, C. Dullemond, [19] 

G. Rupp and J. E. Ribeiro, Z. Phys. C 30, 615 (1986). 
[2] T. Kunihiro, Prog. Theor. Phys. Suppl. 168, 156 (2007). [20 

[3] K. Hagiwara et al. [Particle Data Group], Phys. Rev. D [21 

66, 010001 (2002). 
[4] S. Eidelman et al. [Particle Data Group], Phys. Lett. B 

592, 1 (2004). [22 

[5] I. Caprini, G. Colangelo and H. Leutwyler, Phys. Rev. 

Lett. 96, 132001 (2006). [23 

[6] E. M. Aitala et al. [E791 Collaboration], Phys. Rev. Lett. 

86, 770 (2001). 
[7] M. Ishida, S. Ishida, T. Komada and S. I. Matsumoto, [24 

Phys. Lett. B 518, 47 (2001). 
[8] D. V. Bugg, Phys. Lett. B 572, 1 (2003) [Erratum-ibid. [25 

B 595, 556 (2004)]. 
[9] F. E. Close and N. A. Tornqvist, J. Phys. G 28, R249 [26 
(2002). 

[10] R. L. Jaffe, Phys. Rev. D 15, 267 (1977); D. Black, [27 

A. H. Fariborz, F. Sannino and J. Schechter, Phys. Rev. 
D 59, 074026 (1999). [28 

[11] G. 't Hooft, G. Isidori, L. Maiani, A. D. Polosa and V. Ri- 

quer, Phys. Lett. B 662, 424 (2008). [29 

[12] T. Hatsuda and T. Kunihiro, Prog. Theor. Phys. 74, 765 

(1985). [30 

[13] C. E. Detar and J. B. Kogut, Phys. Rev. D 36, 2828 [31 

(1987); M. G. Alford and R. L. Jaffe, Nucl. Phys. B 578, 

367 (2000); W. J. Lee and D. Weingarten, Phys. Rev. D [32 

61, 014015 (2000); C. McNeile and C. Michael [UKQCD [33 

Collaboration], Phys. Rev. D 63, 114503 (2001). 

[14] W. A. Bardeen, E. Eichten and H. Thacker, Phys. Rev. [34 

D 69, 054502 (2004). [35 

[15] S. Prelovsek, C. Dawson, T. Izubuchi, K. Orginos and 

A. Soni, Phys. Rev. D 70, 094503 (2004). [36 

[16] H. Suganuma, K. Tsumura, N. Ishii and F. Okiharu, PoS 

LAT2005, 070 (2006). 
[17] N. Mathur et al., Phys. Rev. D 76, 114505 (2007). 
[18] T. Burch, C. Gattringer, L. Y. Glozman, C. Hagen, [37 
C. B. Lang and A. Schafer, Phys. Rev. D 73, 094505 
(2006). [38 



C. McNeile and C. Michael [UKQCD Collaboration], 
Phys. Rev. D 74, 014508 (2006). 

K. Hashimoto et al. [RBC Collaboration], in preparation. 
N. Ishii, T. Doi, H. Iida, M. Oka, F. Okiharu, H. Sug- 
anuma and K. Tsumura, AIP Conf. Proc. 842, 492 
(2006). 

H. Wada, T. Kunihiro, S. Muroya, A. Nakamura, C. Non- 
aka and M. Sekiguchi, Phys. Lett. B 652, 250 (2007). 
T. Kunihiro, S. Muroya, A. Nakamura, C. Nonaka, 
M. Sekiguchi and H. Wada [SCALAR Collaboration], 
Phys. Rev. D 70, 034504 (2004). 

M. A. Shifman, A. I. Vainshtein and V. I. Zakharov, Nucl. 
Phys. B 147, 385 (1979); ibid. 448. 

L. J. Reinders, S. Yazaki and H. R. Rubinstein, Nucl. 
Phys. B 196, 125 (1982). 

T. V. Brito, F. S. Navarra, M. Nielsen and M. E. Bracco, 
Phys. Lett. B 608, 69 (2005). 

H. J. Lee and N. I. Kochelev, Phys. Lett. B 642, 358 
(2006). 

R. D. Matheus, F. S. Navarra, M. Nielsen and R. Ro- 
drigues da Silva, Phys. Rev. D 76, 056005 (2007). 
H. X. Chen, A. Hosaka and S. L. Zhu, Phys. Rev. D 76, 
094025 (2007). 

T. Kojo and D. Jido, Phys. Rev. D 78, 114005 (2008). 
V. L. Chernyak and A. R. Zhitnitsky, Nucl. Phys. B 201, 
492 (1982) [Erratum-ibid. B 214, 547 (1983)]. 
S. H. Lee and H. C. Kim, Nucl. Phys. A 612, 418 (1997). 
H. C. Kim and S. H. Lee, Z. Phys. A 357, 425 (1997); 
arXiv:nucl-th/9610013 

D. B. Leinweber, Ann. Phys. (NY.) 254, 328 (1997). 
V. Dmitrasinovic, Phys. Rev. C 53, 1383 (1996); Nucl. 
Phys. A 686, 379 (2001). 

K. Naito, M. Oka, M. Takizawa and T. Umekawa, Prog. 

Theor. Phys. 109, 969 (2003). T. Umekawa, K. Naito, 

M. Oka and M. Takizawa, Phys. Rev. C 70, 055205 

(2004). 

V. Elias, A. H. Fariborz, F. Shi and T. G. Steele, Nucl. 

Phys. A 633, 279 (1998). 

E. V. Shuryak, Rev. Mod. Phys. 65, 1 (1993). 



[39] E. M. Aitala et al. [E791 Collaboration], Phys. Rev. Lett. [40] B. Grinstein, R. Jora and A. D. Polosa, Phys. Lett. B 
89, 121801 (2002); M. Ablikim et al. [BES Collabora- 671, 440 (2009). 

tion], Phys. Lett. B 633, 681 (2006). 



